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I.  INTRODt ICT ION 


2. 1  Hirst  type  heal  equation 


Almost  all  the  theories  about  the  phenomenon 
of  El  Nino  and  Southern  Oscillation  (F.NSO)  as¬ 
sume  a  positive-only  correlation  in  the  ocean  heat 
equation  between  the  time  rate  change  of  sea  sur¬ 
face  temperature  (SST)  anomaly  (flTjdt)  and  the 
upper  ocean  thickness  fluctuation  (/r).  It  is  very 
important  to  justify  the  validity  of  this  positive-only 
correlation  before  developing  complicated  coupled 
ocean  atmosphere  models. 

Based  on  classical  K  rnus  Turner  type  ocean 
mixed  layer  (OML)  theory,  it  is  found  that  (lie 
positive-only  correlation  between  and  h  is 

due  to  overlook  of  an  ultimately  important  feature 
in  the  OKI I.:  iron-upward  turbulent  heal  this  at  the 
OMI  base  when  the  OML  temperature  is  higher 
than  the  temperature  below  the  OML,  as  it  usually 
is  in  the  tropics.  With  consideration  of  this  feature, 
both  positive  and  negative  correlations  between  the 
time  rate  change  of  SST  anomaly  and  OML  depth 
fluctuation  have  been  found  Furthermore,  if  initial 
OML  dcpllr  is  greater  than  the  Monin-Ohukhov 
length  scale  (Z).  the  OML  quickly  shallows  lo  the 
Monin  Oluikhov  length  scale,  and  the  negative 

correlation  exists;  however,  if  initial  OML  depth  is 
smaller  than  tire  Monin  Obukhov  length  scale,  lire 
positive  correlation  appears.  Shift  of  positive  only 
correlation  lo  positive,  negative  correlation  will  lead 
to  a  new  FNSO  theory. 


An  oceanic  heat  equation  oRen  appeared  in 
simple  coupled  ocean  atmosphere  models  (eg. 
Hirst.  1986;  Hirst  and  Lau.  1990)  was  First  pre 
sented  by  Hirst  (1986) 

dT.  _  .  A  A 

-^-  =  -Tji  +  K,h-9tTt  (I) 

where  7*„  h,  u  arc  perturbations  c*F  SST.  tipper 
ocean  thickness,  and  zonal  currents:  T,  is  the  mean 
zonal  temperature  gradient;  »,  is  thermal  dissi 
pation  coefficient;  and  AT,(  >  0)  is  SS  I  thickness 
coupling  coefficient  I  lie  positive  value  of  R,  . 
caused  by  the  only  appearance  of  one  of  the  two 
cases  of  the  OMI  (i.e.,  the  OMI  thickness  less  than 
the  Monin-Ohukhov  length  scale),  guarantees  the 
positive-only  correlation  between  DTjdl  and  h 

As  pointed  by  Xic  cl  al  (1989).  tire  SSd 
equation  of  Anderson  and  McCreary  (|9X5)s  cou¬ 
pled  air-ocean  m«»dcl  also  belongs  to  this  category. 


2.2  Zcbiak-Canc  type  heat  equation 

In  several  coupled  numerical  models  (e.g.. 
Zebiak  and  Cane  1987;  Baltisli.  1988;  Ratlisti  and 
Hirst.  1989;  Neelin  1991).  the  following  SST 
anomaly  equation  is  often  used: 


2  SSI  EQUATION  IN  CURRFNT  FNSO  1111 
OKIES 


Cy—  A(*r 


4  w) 


In  the  current  ENSO  theories  containing 
ocean  thermodynamics,  almost  all  imply  a 
positive  only  correlation  between  lime  rale  change 
of  SSI  anomaly  and  OML  depth  fluctuation  I  cl 
us  investigate  the  two  typical  (Hirst  type  and 
Zebiak  Cane  type)  heat  equations. 


-  {Ali?  4  vv)  -  AI»T)|  xf?-  «,f,  (2) 

where  CT  is  the  horizontal  temperature  advection; 
//,  is  the  mean  upper  layer  thickness;  ie,  u-  are  mean 
and  fluctuation  of  vertical  velocity;  and  T.„h  is  the 
subsurface  temperature,  which  is  determined  em¬ 
pirically  by 


52 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

1992 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1992  to  00-00-1992 

4.  TITLE  AND  SUBTITLE 

An  Important  Ocean  Feature  Overlook  in  Current  El  Nino-Southern 
Oscillation  Theories 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Department  of 

Oceanography , Monterey  ,C  A  ,93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Third  Conference  on  Global  Change,  American  Meteorological  Society,  52-59 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

1 

1 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 
ABSTRACT 

Same  as 
Report  (SAR) 


18.  NUMBER 
OF  PAGES 

8 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Tu,b  =  ^il  lanh  l*|(*  »  A)|  — h  >  0 

^isuh  =  7j|  •an*1  l^(A  —  /»)]—  lnnli(/>i/»)|.  /;  <  0 
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where  the  coefficients  m  Zebiak  C'nne  model  take 
the  following  values:  r  =  0.7  .  T,  =  2 XT. 

F,  =  -  40' C  .  b,  =  (XOm!  '  .  and  ki—(33mjr*. 
Rattisti  ( I OXS )  ,md  Ratlisli  and  Hirst  (I9S9)  uscil  a 
little  dilTercnl  form  for  7,.h: 

7Mlh  =  </(/!  )/i.  ,i(Ti)>0  (3/0 

Recently.  Ncelin  (1991)  adopted  a  linear  form  for 
7^: 

7jqh  =  ®o  +  *i^  (VI 

where  »„  and  a,  are  two  empirical  constants.  I  he 
function  A(.v)  in  (2)  is  defined  by 

A(.v)-0.  ,ts;<1 

A(.r)  —  x.  x  >  0  (4) 

Utilization  of  any  one  of  (3a.b.c|  with  (2)  leads  to 
a  positive-only  correlation  hetween  0 TJSt  and  h. 
i  c..  the  increase  (or  decrease)  of  h  causes  the  in¬ 
crease  (or  decrease)  of  STJSt.  Is  this  positive  only 
correlation  correct?  I  he  answer  is  no.  Recently, 
Chu  (I99la.b,  c)  pointed  nut  the  possibility  of  ev 
ixtcivce  of  negative  correlation  between  STJSt  and 
h  base  on  the  surface  wind  condition  In  this  note, 
a  further  investigation  of  thermodynamical  and 
dynamical  structure  of  the  OMI  explores  tlie 
oceanic  condition  for  the  negative  correlation. 

3.  OCEAN  MIXED  LAYER  MODI  I 

Roth  almosphcric  and  oceanic  planetary 
boundary  layers  are  ultimately  important  (or 
ocean-atmosphere  interaction  since  fluxes  of  mo¬ 
mentum  and  buoyancy  across  the  air  ocean  inter¬ 
face  are  major  driving  forces  (or  dampening 
factors)  fot  the  two  fluids.  Therefore,  before  do 
signing  any  I.NSO  model  (either  theoretical  or  mi 
mcrical),  whether  the  boundary  layer  physics  is 
reasonable  or  not  should  he  first  checked  up 

Kraus  and  Turner  (l%7)  developed  a  OMI. 
model  by  considering  the  layer  to  lie  vertically  ho 
mogeneous:  heat  inputs  to  the  layer  ami  mass  cn 
trained  at  live  bottom  of  the  layer  are 
instantaneously  mixed  uniformly  through  the  layer. 


Turbulence  in  the  upper  ocean  is  strong  to  main¬ 
tain  a  vertically  uniform  temperature  down  to  a 
certain  depth,  h  The  OMI  temperature  (T.)  and 
dcplh  (/>)  arc  important  variables  predicted  by 
OMI  models.  Rclow  the  depth  h.  there  is  a  rapid 
decrease  in  temperature  across  what  is  called  the 
entrainment  zone,  live  thickness  of  entrainment 
zone  r  is  assumed  infinitesimally  small  Turbulence 
is  assumed  to  he  zero  below  the  entrainment  zone, 
where  the  vertical  velocity  tr  is  determined  by  the 
dynamics  of  the  ocean  interior.  The  model  OMI 
is  depicted  in  fig. I.  Vertically  uniform  temperature 
in  the  OMI  (therefore,  no  vertical  advcction  in  the 
OMI.)  and  no  turbulence  in  the  deeper  layer  below 
the  entrainment  zone  lead  to  three  different  forms 
of  heat  equation  for  three  parts  of  the  ocean 

ST.  a  -  n  c 

(wT I  =  £  -h  .  <„i 

St  oz  ilz 


i n 

St 


+  if  —  t  ^(w'T) 

07  02 


It 


(5 b) 


ST  ST  SF 

St  Sz  Sz 


z  <  —h  —  t  (V) 


where  w'T  is  vertical  turbulent  heat  flux;  is  (he 
characteristic  value  of  sen  water  density;  cpw  is  spe 
cific  heat  of  sea  water  under  constant  pressure,  and 
T  is  downward  flux  of  solar  radiation  in  the  ocean 
computed  by 


is 


F„  = 


Ip 

I’m  ifpw 


(6) 


where  /„ ,  y  are  the  solar  radiation  incident  on  the 
ocean  surface  and  the  average  extinction  coefficient. 
Vertical  integration  or  heal  equation  (5)  from  the 
•Kean  surface  to  the  OMI  base  leads  to 


L= <w7r)_*  -  (wrv,  +  f o  -  f  * 


(7) 


At  the  ocean  surface,  the  turtvulent  heat  flux  must 
equal  the  net  heat  transfer  through  the  air-ocean 
interface,  i.e.. 


(w‘’7*V,  = 


Rh  +  IT  +  i.e 


Pu*  tfpw 


(XI 


where  Rt  is  the  net  heat  loss  by  longwave  radiation 
from  the  sea  surface;  and  II,  and  I.E  are  the  upward 
sensible  and  latent  heat  fluxes  across  live  air  ocean 
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interface.  Vertically  integrated  turbulent  kinetic 
energy  (TK It)  equation  is  given  by  (Kraus  and 
Turner,  1 967) 


Elimination  of  dTjdt  from  (5)  and  (II)  anti  utili 
zation  of  (9a)  leads  to  a  relation  for  computation 
of  turbulent  heal  flux  at  the  OML  base 


w'T<1z  =  -  (G  />) 


<r>r/) 


where  D  denotes  the  dissipation  within  the  OMI 
and  G  the  TKE  input  from  the  wind.  The  simplest 
parameterization  for  (G  —  D)  is  (Denman,  1973) 


0-/)=  -  > n 


(9/,) 


("•-n  h  =  -(n'V\}- 


2 (G  -  D) 
It 


+  2[  ( I  4-  e-1^r)  — ( I  —  e^1)]  (12) 

If  the  penetration  depth  of  solar  radiation  is  much 
smaller  than  the  OML  thickness,  i.c.,  j'/nH,  the 
exponent  r  ,h  is  negligible  against  I,  and  therefore 
Eq.(l2)  can  be  simplified  as 


where  V„  is  the  horizontal  wind  at  It)  m  height: 
is  the  characteristic  value  of  surface  air  density:  C„ 
is  the  drag  coefficient;  g  is  gravity;  a  is  thermal  ex¬ 
pansion  coefficient  of  sea  water;  and  m  is  propor¬ 
tionality.  Integration  of  the  heat  equation  (.5)  with 
respect  to  z  from  a  depth  -  inside  the  OML  ( 
z  >  —It)  to  the  ocean  surface  leads  to 

dr. 

+  tw’7'% - (i/n*  -  ru -  ro  i mi 

Further  integration  of  (10)  with  respect  to  z  from 
the  OML  base  to  the  ocean  surface  gives  rise  to 


Uv'Vlj,  =  —  2[  \  (i*'T)(,  -T  _  ru(  |  ±  >| 

(13) 

The  turbulent  heal  flux  at  the  OML  base  is  either 
downward  or  zero,  never  upward,  i.c., 

0^n_,,<;0  (M) 

due  to  higher  OML  temperature  (T.)  versus  deeper 
layer  temperature  (Fig.l).  With  consideration  of 
this  constrain,  a  better  form  is  recommended  to 
replace  ( 13) 


*iiZi 

2  !U 


ft) 

Vrdz  -  h(V'T\  +  IK  r()  + 


(II) 


(tv'V)  „  =  - 


C?u 


P»o<) 


A( 


r  -  h 


I'"1 


(15) 


where  Qu  is  the  net  downward  heat  flux  plus  radi¬ 
ation  at  the  ocean  surface; 


Qn  -  Fiv-n'/nJ  l(t  ~  ( l'\)\  (I  &) 

which  is  generally  positive  in  the  (topics  when  the 
low  frequency  mode  is  concerned;  and  f  is  Monin 
Obukhov  length  scale  defined  by 


2(0-1)  -I  /’()/ )' ) 


(17) 


Integrating  the  heat  equation  (5b)  across  the 
entrainment  zone  from  z  —  —  In  ;=  -  It  and 

taking  the  limit  as  r.i  >f),  a  jump  condition  has  been 
obtained  (Denman.  1973) 


(w'  -rU  =  -A(iv  /1  +  ^  )(7C  -  7 .,,)  (18) 


where  w  T_h  arc  vertical  velocity  and  temperature 
at  (he  base  of  the  entrainment  zone,  respectively. 


Fig  I.  Ocean  mixed  layer  model. 
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In  fact. 


,  dh 

7,+  ^ 


is  called  the  en  trainmen  I  velocity  which  were  used 
in  some  I.NSO  models  containing  the  OML  (e  g 
Anderson  and  McCreary,  1985;  Hirst,  |98(*;  Xie  et 
nl,  1989).  Since  (lie  OMI  temperature  is  higher 
than  the  temperature  below  the  OMI  ..  Ihc  function 
A(t)  alv>  guarantees  non  upward  turbulence  heal 
flu*  at  the  OMI  base.  If  w  ,  +  i'lltJPi  >  0.  there  is 
entrainment  mixing  (causing  downward  turbulent 
heat  flux)  at  the  base  of  the  entrainment  /one.  II 
w  +  dftfflt  £  0,  (here  is  no  entrainment  mixing 
(causing  zero  turbulent  heal  flux)  at  the  Ixasc  of  tl*c 
entrainment  zone.  Fqs.(7),  (15),  (17).  and  (18)  are 
basic  equations  for  Kraus  I  urner  type  OMI.  mod 
els  Eq.(IS)  is  used  to  compute  turbulent  heat  flux 
(w'V)^,  Fq.(7)  is  used  In  determine  the  OMI 
temperature  T„  and  Fq.(IS)  is  used  to  calculate 
lime  rate  change  of  OML  depth  ilhfdt  in 
entrainment  mixing  case,  e  g.. 


dh  WrU 


l/r  <  /■>  (19) 


From  Wyrtki's  ( 1981)  estimation.  Ihc  mean  upwell 
ing  velocity  at  the  base  of  the  Kkman  layer  ( 
7  =  —  ft,  )  Ls  nearly  I  mftkty.  Theoretically.  tl»c 
F.kmnn  depth  at  the  equator  is  infinity.  It 
<i,~200m  is  thought  to  be  a  reasonable  estimation 
of  Kkman  layer  thickness,  the  mean  vertical  velocity 
at  (lie  mean  nmed  laser  depth  can  be  roughly 
evaluated  by 

-r—  iv  /jiO.  I  m/itay 

'’/i 

where  Ihc  mean  mixed  layer  depth,  taken  as  24  lit. 
is  from  ihc  data  collected  during  Hawaii  lo-Tahiti 
Shuttle  E  xperiment  (.Schneider  anil  Muller.  1990). 

Fig. 2  indicates  the  OMI  deepening  rate 
caused  by  the  turbulent  beat  llux  at  the  OMI  base 
Tor  fixe  different  values  of  the  temperature  jump  at 
the  OMI  base  (7‘  —  T  ,1  For  h  not  veil  close  to 
f .  in  the  righthand  side  of  (19)  the  second  term  is 
much  larger  than  the  lirsl  term,  which  means  lhal 
Ihc  change  of  OMI  depth  is  more  likely  by  Ihc 
strength  of  turbulent  heal  llux  al  Ihc  OMI  base 
rather  than  by  ihc  Kkman  pumping  Eq.(9l>)  indi¬ 
cates  that  an  increase  nl  surface  xvind  speed  trails 
to  an  increase  of  (•  —  />.  which  in  turn  augments 
the  OMI  deepening  rale  through  the  increase  of 


the  downward  turbulent  heal  flux  at  the  OMI  base, 
—  (w'T)  as  It  <  f  | sai 1 3)| ;  or  increases  the  OMI 
depth  as  h»f  |.wr(l7)|.  In  othci  u’nrds.  the  in 
crease  of  surface  wind  speed  generally  increases  the 
OMI  depth 


fig. 2.  Dependence  of  OML  deepening  rale  (for 
fi  <  O  caused  by  tire  turbulent  heat  flux  al  the 
OML  base.  Ur"D_J{T,  —  T on  /i/f  for  different 
values  of  (7,-  F,)  :  (a)  IT.  (b)  2*C.  (c)  VC  (d) 
4T.  and  (e)  5T. 


4  AN  IMPORTANT  FEATURE  OVIK 
LOOKED  IN  TIIF  ENSO  Til  KORIES 


Another  form  of  OMI.  beat  equation  can  l>c 
obtained  by  the  substitution  of  (15)  into  (7): 


I 


fli 

P» 


(  -It 
h 


)|(20| 


Here.  is  the  characteristic  OMI  thickness.  Hie 
time  rale  change  of  SSI  is  determined  by  ihc 
Monin Obukhov  length  scale  f  and  OMI  depth  It 


An  ultimately  important  feature  has  been 
overlooked  in  tire  current  ENSO  theories  is  the  ex¬ 
istence  of  two  different  analytical  forms  of  the 
OMI  heal  equation  depending  on  Ihc  initial  OMI 
drplh  If  lire  inilial  OMI  depth  is  greater  Ilian  the 
Monin  Obukhov  length  scale  (A>f).  the  computed 
value  of  the  turbulent  heat  flux  al  the  OML  base 
(iiT]4,  which  is  only  determined  by  the  ocean 
surface  wind  forcing  (fi  />).  OMI  depth  (fi),  and 
the  surface  ncl  heal  flux  |trr(l1)|.  will  lie  positive, 
which  is  against  tire  physical  law  lhal  there  is  no 
upward  turbulent  heat  flux  at  the  OML  (rase  xvtren 
T,  >  T  t  Moreover,  the  function  A((/’  h)/h\  -  0. 

which  makes 
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(h  -r)  (21//) 


(^r)/(  =  n 

and  (urns  Ihc  SST  equation  (20)  into 

r  Q)  |  I  _*%  _  jf 

PsvOt:,m^  dt  >'■ 


which  shows  that  the  time  rate  change  of  SSI  ( 
PTjtU)  monotonically  decreases  with  the  increase 
of  h  (Fig.3).  It  is  easy  to  find  in  this  case  that  the 
increase  of  the  OML  depth  h  causes  Ihc  decrease 
of  the  time  rate  change  of  SSI.  which  implies  a 
negative  correlation  between  Ihc  time  rate  change 
of  SST  anomaly  and  OML  deplh  fluctuation,  which 
was  overlooked  in  the  current  LIS'SO  models. 


1  'g-4.  Time  rate  change  of  SS  I  versus  OML  thick 
ness  for  h  >  f. 


fhg-  L  I  inie  rate  change  of  SS  I 
ness  for  h  <  r. 


versus  DM  I 


I  hick - 


As  the  OML  depth  is  smaller  than  ihc 
Monin  Obukhov  length  scale  (//</).  i.e.,  Ihc  lur 
bulenl  heal  flux  at  I  lie  OML  base  is  negative 

wr  u  <  o 


and  the  OML  heat  equation  (20)  becomes 


a, 


!)t 


=  {h  <0(2l/») 


which  shows  that  the  lime  rate  change  of  SST  ( 
flTJrU)  monotonically  increases  with  Ihc  increase 
of  h  (Fig.4).  This  leads  to  a  positive  correlation  he 
tween  the  time  rate  change  oT  SST  anomaly  and  Ihe 
OML  depth  fluctuation,  which  has  been  well  re- 
presented  in  the  current  FNSO  theories. 


5.  GENERATION  OF  TWO  TYPES  OF  UN- 
STABLE  MODES 

Since  the  OML  heat  equation  has  two  differ¬ 
ent  forms  (21a)  and  (21b)  due  to  initial  OML  depth 
versus  the  Monin-Obukhov  length  scale,  it  is  antic¬ 
ipated  that  the  features  of  unstable  modes  gener¬ 
ated  by  coupled  air- ocean  system  should  he 
different. 

The  tropical  subcloud  layer  is  usually  domi¬ 
nated  by  mean  easterlies.  The  low-level  atmo¬ 
spheric  convergence  and  convection  ('Convection 
A")  arc  generally  associated  with  the  warm  water 
(Gill,  1980).  In  ihc  cast  (or  west)  ol  Ihe  warm  water 
region,  the  surface  wind  speed  is  enhanced  (or  re¬ 
duced).  and  therefore,  (he  OMI  deplh  h  is  in¬ 
creased  for  decreased).  There  are  two  main  cases 
which  arc  illustrated  as  follows. 

5.1  Small  initial  OML  thickness  (/r  <  f  ) 

In  this  case,  the  SSI  equation  is  given  by 
(21b),  showing  a  positive  correlation  between 
PTJdl  and  the  OML  depth  (Fig. 5).  Therefore,  a 
positive  OML  fluctuation  (//  >  0)  leads  to  a  positive 
time  rate  change  of  SS  I  anomaly  >  0).  In 

the  east  (nr  west)  of  Convection  A.  the  OML  depth 
fluctuation  is  positive  (or  negative),  which  in  turn 
gives  rise  to  a  positive  (or  negative)  time  rale 
change  of  SS  I  anomaly  (Fig. 5)  I  lui.s,  positive  SS  I 
anomaly  and  associated  unstable  convective  dis¬ 
turbances  (denoted  by  'Convection  B')  will  appear 
in  the  east  of  the  initial  warm  water  region.  Ap 
plying  the  same  argument  to  Convection  B,  positive 


56 


I 


SSI  anomaly  and  associated  unstable  convective 
disturbances  (denoted  by  "Convection  C")  will  ap¬ 
pear  in  the  east  or  Convection  B  (I  ig,f>)  due  to  Hie 
rnrtber  increase  of  OMI  depth  caused  by  (he  cm 
lianccd  surface  easterlies  (sum  of  l lie  background 
easterlies,  easterlies  associated  with  Convections  A 
and  B).  Unstable  convective  disturbances  arc  un¬ 
likely  to  appear  between  Convection  A  and  Con 
vection  B  due  to  the  opposite  directions  of  the 
surface  winds  associated  with  A  and  B  Therefore, 
as  li  <  f ,  the  unstable  mode  is  generated  in  Hie  east 
ol  the  initial  warm  water  region  and  propagates 
eastward. 

For  a  typical  La  Nina  condition,  atmospheric 
deep  convection  (Convection  A)  develops  over  Hie 
western  pacific  warm  water.  A  shallow  and  cool 
OMI  water  occupies  the  eastern  Pacific.  The  OMI 
thickness  is  generally  satisfied  the  condition  h  <  / . 
The  warm  SS  I  disturbances  and  associated  atmo¬ 
spheric  convection  arc  generated  in  the  east  oT 
Convection  A  and  propagate  eastward,  which  leads 
to  a  eastward  propagation  of  t lie  surface  westerlies 
(westerlies  associated  with  Convections  A,  B,  and 
C).  At  (be  same  lime,  the  OMI.  thickness  h  in¬ 
creases.  The  whole  process  features  Ihc  transition 
from  La  Nina  to  El  Nino. 


tiTjUt  -  h  <  II 

_jQ _ 

OML  thickness 


H  »  n 


TJ - 

OMI-  thickness 
increases 


Fig  5.  Positive  SST  anomaly  and  associated 
conveclive  disturbances  (Convection  B)  generated 
in  the  east  of  Convection  A  (h  <  /) 


OML  tlurlnrii  CJ- 

decreases  OML  thickness 


increases 


Fig. 6.  Eastward  propagation  of  positive  SSI 
anomaly  and  associated  convective  disturbances  ( 
h  <{).  ' 


5.2  Large  initial  OML  thickness  (/)  >  f) 

In  this  case,  the  OML  quickly  shallows  to  the 
Monin  Obukhov  length-scale.  The  OML  shallowing 
lime  scale  is  estimated  by  (Garwood,  1 077) 


where  £  is  the  mean  TKF.  in  the  OML.  If 
<  E  >,tI  ~  10  rm/.t  and  li  ~  100 in,  the  OMI.  shal¬ 
lowing  time  scale  (t)  is  nearly  1000  seconds,  which 
is  so  short  comparing  to  the  ENSO  lime  scale. 
Thus,  it  is  reasonable  to  assume  that  the  OMI.  in¬ 
stantly  shallows  to  the  Monin-Obukbov  length  scale 
/  when  the  initial  OML  depth  is  greater  than  f. 
The  SSI  equation  is  given  by  (21a),  showing  a 
negative  correlation  between  ctTjfit  and  the  OML 
depth  (I ;ig.2),  which  means  that  a  negative  OMI 
fluctuation  (//  <  0)  leads  to  a  positive  time  rate 
change  of  SST  anomaly  (<77*, /tV  >  0).  In  the  west 
(or  east)  of  initial  Convection  D,  the  OMI  depth 
fluctuation  is  negative  (or  positive),  which  in  turn 
gives  rise  to  a  positive  (or  negative)  time  rate 
change  of  SST  anomaly  (Fig. 7)  I  litis,  positive  SST 
anomaly  and  associated  unstable  convective  dis¬ 
turbances  (denoted  by  "Convection  E  )  will  appear 
in  the  west  of  Convection  D.  Applying  (lie  same 
argument  to  Convection  E,  positive  SSI  anomaly 
and  associated  unstable  convective  disturbances 
(denoted  by  "Convection  F")  will  appear  in  the  west 
of  Convection  F  (Fig.X)  due  to  Ihc  further  decrease 
of  OMI  depth  caused  by  the  reduced  surface 
easterlies  (the  background  easterlies  cnunlerhal 
anced  by  westerlies  associated  with  Convections  F> 
and  E).  Unstable  convective  disturbances  arc  un 
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likely  to  appear  between  Convection  I)  and  Con 
vection  !■  due  to  the  opposite  directions  of  die 
surface  winds  associated  with  I)  amt  FI  T  herefore, 
as  h  >  f .  the  unstable  mode  is  generated  in  the  west 
of  the  initial  warm  water  region  ami  propagates 
westward. 

l  or  a  typical  El  Nino  condition,  a  deep  ami 
warm  OML  water  occupies  the  eastern  Pacific, 
where  llie  atmospheric  deep  Convert  inn  (Con¬ 
vection  I))  develops  The  large  OMI  thickness  is 
quickly  adjusted  to  the  Monin  Obukhov  length 
scale,  h  =  f  the  warm  SSI  disturbances  ami  asso¬ 
ciated  atmospheric  convection  are  generated  in  llie 
west  of  ( onvcclion  D  and  propagate  westward, 
which  leads  |o  a  westward  propagation  of  tin*  sur 
face  westerlies  (westerlies  associated  with  Con 
vections  I).  E.  and  F).  At  the  same  lime,  the  OMI 
thickness  h  decreases.  The  whole  process  features 
the  transition  from  El  Nino  to  l  a  Nina 


Fig, 7.  Positive  SST  anomaly  and  associated 
convective  disturbances  (Convection  F.|  generated 
in  llie  west  of  Convection  D  (k  =  f). 


Fig.X  Westward  propagation  of  positive  SST 
anomaly  ami  associated  convective  disturbances  ( 


6.  SUMMARY 

(1)  Both  eastward  and  westward  propagating 
unstable  modes  described  here  is  produced  purely 
by  the  SST-OML  thickness  feedback.  Positive  (or 

A  A 

negative)  correlation  between  P'I'JPt  and  h  will 
generate  eastward  (or  westward)  propagating  un 
stable  mode. 

(2)  Shift  from  positive  only  correlation  to 
posit i ve/negative  correlation  between  PI, /Pi  and  It 
leads  to  a  discovery  of  a  new  mechanism  of  La 
Nina  and  El  Nino  transition  based  on  the  magni¬ 
tude  of  OML  thickness  (/«)  relative  to  the  Monin- 
Obukhov  length-scale  (/)  It  also  indicates  the 
existence  of  an  OML  switcher  for  this  transition. 

(3)  The  vertical  advert  ion  Icon  can  not  ex 
plicitly  appear  in  the  OML  heal  equation  because 

it  only  afTects  the  turbulent  heal  Hus  at  the  OML 
base  through  lire  jump  condition  (IX).  As  h  <  f.  the 
vertical  advert  inn  changes  tire  OML  temperature 
through  the  change  of  OML  thickness  As  h'Z.f. 
the  OML  quickly  shallows  to  the  Monin  Obukhov 
length  scale  and  llie  turbulent  heal  flux  at  lire  OML 
base  heconrcs  zero.  Therefore,  in  this  case  lire  ver¬ 
tical  advection  no  longer  affects  lire  OMI 


ACKNOWLEDGE  MEN  I  S 


This  work  was  funded  by  the  OITice  of  Naval 
Research  and  the  Naval  Postgraduate  School.  The 
author  thanks  ILL.  Ilaney  ami  R.W,  Garwood.  Jr. 
at  the  Naval  Postgraduate  School  for  invaluable 
discussion. 


58 


REFERENCES 

Anderson.  0.1. .T..  and  J.P  McCreary,  I9RV 
Slowly  propagating  disturbances  in  a  coupled 
ocean  atmosphere  model.  J.  A  linos.  Sci.,  42. 
615  629 

Ralltsti.  D.S .  1988:  I  he  dynamics  and 

thermodynamics  of  a  warming  event  in  a  coupled 
tropical  a  I  mosphe  re-«  icean  model  J.  A  linos.  Sci.« 
45,  >889  >9 |9. 

Raltisli.  D.S.,  and  A  C  Hirst.  1980  Inlerannual 
variability  in  tire  tropical  atmosphere  ocean  system 
influence  of  the  hasic  stale,  nonlinearity,  and  ocean 
geometry  J.  Atmos.  Sci..  46.  1687-1  12 
Chu,  PC..  1991a  Generation  of  low  frequency  tin 
stable  modes  in  a  coupled  cqwaloiial  troposphere 
and  ocean  limed  layer  model  J.  Atmos.  Sci_. 
submitted 

Chu.  PC..  I99|b:  Air-ocean  surface  heat  exchange 
model  and  low  frequency  unstable  modes  in  at 
ntosphcrc  and  ocean.  Proceedings  on  Fifth  Con¬ 
ference  on  (  lunate  Variation,  Am  Meleond.  Sot 
in  press. 

Chu,  P  C..  |99lc:  Is  ocean  Ihermndynnmics  correct 
in  current  El  Nino  and  Southern  Oscillation  llico- 
ries?  I  Gaoohvs.  Res.,  submilted. 

Dennutn,  K  i ...  1973:  A  time-dependent  model  bl 
the  upper  ocean  J.  Phys.  Occnnogr.,  3,  171  184. 
Clill,  A.E.,  1980:  Some  simple  solutions  for  heat- 
induced  tropical  circulation.  Quart.  J,  Roy. 
Meteor.  Soc.,  106,  447-462. 

Hirst.  A.C  .  1986:  Unstable  and  damped  equatorial 
modes  in  stmplc  coupled  atmosphere-ocean  models. 
J.  Atmos.  Sci.,  43.  606  630. 

Hirst.  A.C  .  ami  K  M  Lait.  1990:  Inlraseasonal  and 
intra-annual  oscillations  in  coupled  ocean- 
atmosphere  models.  J.  Clim.,  3,  713-723 
Kraus,  E.B.  and  J.S.  Turner.  1967:  A  onc- 
drmensional  model  of  the  seasonal  Ihermocline  II 
The  general  theory  and  its  consequences.  Icllus, 
12,  98  106. 

Ncelin.  I  I)..  I99|:  The  slow  sea  surface  temper¬ 
ature  mode  and  the  fast-wave  limit:  analytic  theory 
for  tropical  inlerannual  oscillation  and  experiments 
in  a  hvbrid  coupled  model.  J.  Atmos.  Sci.,  48, 
584  606 

Schneider.  N.  and  P  Muller,  1990;  "Hie  meridional 
and  seasonal  structures  of  the  mixed  layer  depth 
and  its  diurnal  amplitude  obscivcd  during  Hawaii 
In-Tahiti  shuttle  experiment  J.  I'hys  Occanogr  , 
20,  1395  1404 

VVyrtki,  K..  1981:  An  estimate  of  equatorial  up 
welling  in  Ihe  Pacific.  J.  Phys.  Occnnogr.,  II, 
1205-1214 


Xie,  S.P..  A.  Kubokawa,  and  K  llanawa.  1989: 
Oscillations  with  two  feedback  processes  in  a  cou¬ 
pled  ocean-atmosphere  model  J.  Clim.,  2,  946-964 
7.ebiak.  SC.,  and  M  A  Cane,  1987:  A  model  El 
Nino  Southern  Oscillation.  Mon  Wca.  Rev.,  115, 
2262-  227S. 


59 


